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Abstract 

The classical equations of irrotational water waves have recently been 
reformulated as a system of two equations, one of which is an explicit non- 
local equation for the wave height and for the velocity potential evaluated 
on the free surface. Here we first extend this formalism to n-dimensions, 
n > 2, and then derive rigorously the linear limit of these equations. 
Furthermore, for n = 2, we generalise the relevant forumlation to the case 
of constant vorticity and to the case where the free surface is described 
by a multi-valued function. Also, in the two dimensional case we derive 
a sequence of Hamiltonian systems, hence providing an approximation in 
the asymptotic limit of certain physical small parameters. 
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1 Introduction 



We consider the classical problem in hydrodynamics concerning the propagation 
of surface waves generated by an incompressible fluid with free surface. In the 
case of irrotational flow we consider the problem in n spatial dimensions, whereas 
in the rotational case we confine attention to n — 2. We consider the problem in 
its full generality so that the depth of the fluid may not be constant, and consider 
all solutions, not only those that arise under the travelling wave assumption. 
In addition, we include the effect of surface tension so that capillary-type waves 
are included in our study. 

We denote by Bh the bottom surface: 

B h = {(x,y) : x G R n "\ y = -h + h(x)}, (1.1) 

where ho is constant and h(x) is a real valued function; the notation B^ will 
denote infinite depth. We denote by S v the free surface: 

S v = {(x,y) : x G R"~\ y = r]{x,t)} for t > 0. (1.2) 

We refer to rj as the height of the wave and we assume rj + ho > h for each 
x E R n_1 . The domain of the problem is the region between Bh and S v , which 
is denoted by f2. 



1.1 The Irrotational Case in n Dimensions 

In the irrotational case we introduce the velocity potential 0, where u = V0, 
and then the governing equations become: 

(1.3a) 
(1.3b) 

V.)..V S = )/, nuS, r (1.3c) 

(1.3d) 
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in Q, 





on B h , 


m 


on S v , 


m 


on S v , 



where g is the acceleration due to gravity, V denotes the usual vector gradient, 
Ns is the exterior normal to Bh and Ns is the exterior normal to S n , i.e. 

V = (d x ,d y ), N B = (-d x h,-l), N s = (-d x r),l). (1.4a) 

The right hand side of (jl.3dp is functionally dependent on r\ through: 

P \V 1 + \dxV\ 2 J 

which is a measure of the effect of surface tension (a and p denote the constant 
surface tension and density respectively). Equation (|1.3a|) is a consequence of 
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incompressibility, (|1.3b[) is the Neumann condition satisfied on Bh and (|1.3c[) , 
(jl.3dp are the kinematic, Bernoulli conditions respectively on S v . 

We will reformulate the problem in terms of the functions r\ and q where 
q(x,t) = (j)(x,r],t), i.e. q represents the potential on the free surface S^. It was 
shown by Zakharov in his classical paper (see [T7] ) that in the case n = 2 the pair 
(rj, q) constitutes a canonically conjugate pair in the Hamiltonian formulation of 
this problem. The generalisation of Zakarov's result to n-dimensions is straight- 
forward. Indeed, in the case of constant Bh and zero surface tension equations 
p.3[) admit a Hamiltonian formulation with respect to the Hamilontian 

H = If 5|V</>| 2 dxdy+ / \grfdx (1.5) 
J Jn J 

and the standard symplectic structure. Here dx — dx% A • • • A dx„_i is the 
Lebesgue measure on R 71-1 . 

Applying the chain rule to the expression q(x,t) = <f>(x,r],t) we find the 
following relations: 

d x q = d x (j) + (d y <f>)d x i], (1.6a) 
q t = fa + <t>tVt- (1.6b) 

It is possible to solve for V</> in terms of the q and r/ following a similar calculation 
with that of pQ. Indeed, using (jl.3c[) in equations (jl.6a|) we find the following 
nonsingular set of equations for d x <fi: 

d x q - T]td x r] = (I + d x rj ® d x rf) ■ d x cj> (1.7) 

and then equation p.3cp gives <j) y in terms of q and ~q. In the new variables the 
dynamic boundary condition on S v , i.e. (|1.3d|) . becomes 



In what follows, the notation J dx will be used to denote an integral over R"^ 1 . 



1.2 The Rotational Case in Two Dimensions 

In this section we derive the equations for the free surface problem for a two 
dimensional fluid with constant vorticity. Denoting the velocity of the flow by 
(u,v), the Euler equations for invisid flow are the following equations: 

u t + uu x + vu y = —P Xl (1.9a) 
v t + uv x + vv y = -P y - g, (1.9b) 
u x + v y — 0. (1.9c) 
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Let u> denote the vorticity of the fluid, i.e. 

u = v x -u y . (1-10) 
Eliminating P from (|1.9al) and (|1.9bj) we find: 

^ + (Ud x + Vdy)uj = {). (1.11) 

We restrict attention to the case in which lj — 7 is constant throughout Q, so 
that (|1.11|) is satisfied identically. The domain is simply connected, therefore 
we can introduce a globally defined stream function tp(t, x, y) so that: 

u = ijj y , v = -tp x , (x,y)en. (1.12) 

Replacing in equation (jl.lOp u by ip y and v by — ip x we find Aip — —7. 
Thus, the function defined by ip h = tp + \y 2 l, is harmonic in il. In fact, ip h 
is only defined upto the addition of an arbitrary function of time, so by abuse 
of notation we writeg to mean the [ip ], i.e. the equivalence class of such 
functions that differ only by a function of time. Let tp to be the harmonic 
conjugate of if) h , i.e. 

tp x = ip y + ojy, (1.13a) 
ip y = -tp x . (1.13b) 

The function tp is harmonic throughout J7. 

Adding and subtracting in equations (|1.9a[) and (jl.9bl) the terms vv x and 
uu y respectively, the equations become 

u t + ±d x (u 2 + v 2 )-jv = ~P x , (1.14a) 
vt + k d v( u2 + y2 ) + ^ u = - p v - 9- (L14b) 
Using in equations p,14[) the identities 

Ut = [lpy)t = (<Px)t, V = -ljj x , V t = -(lp x )t = (<Py)t, « = Ipy, 

equations (|1.14p can be integrated to give the following equation: 

Vt + \\VM 2 +ii> + P + gy = a{t), {x,y)eQ, (i.i5) 

where a(t) is some function of time. Since we are still dealing with an equiva- 
lence class of functions tp and tp we can absorb function of time into them, so 
we arrive at 

^ + !lW>| 2 +7^ + .gy = P atm -P, (x,y)eQ, (1.16) 
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where P a tm is the atmospheric pressure above the free surface S v . On the free 
surface S v we have the dynamic boundary condition 

p atm -p=m, m^at^t—) . (i.i7) 

The streamfunction ip is not constant on streamlines in the case of constant 
vorticity, however since we are dealing with an equivalence class of functions we 
may absorb this function of time into ip and use the normalisation ip = on S^. 
Hence (|1.16[) and (| 1 . 1 T[) imply the following nonlinear boundary condition: 

d tV > + ±\Vip\ 2 +gr) = f{v) onS v . (1.18) 

The kinematic condition on and the condition on Bh are the same as in the 
irrotational case. 

In summary, if the vorticity equals the constant 7, then the functions f(x, y, t) 
and T](x, t) satisfy the following boundary value problem: 







Aip 


= 


in CI, 


(1.19a) 


{fx - 


IV, fy) 


•N B 


= 


on Bh, 


(1.19b) 


(fx - 


IV, fy) 


■N s 


= Vt 


on S v , 


(1.19c) 


ft + \fy + \{fx 


- ivf - 


^gy 


= m 


on S v , 


(1.19d) 



where TVg, Ns and /(ry) are defined in (|1.4|) . 

As in the irrotational case, we will formulate the problem in terms of 77 and 
£ where £,(x,t) = ip{x,ri,t), i.e £ represemts the pseudo-potential ip on the free 
surface. Computing £ a using the chain rule (compare with (|1.6aj) ) and invoking 
(|1.19cj) we find the following pair of equations for (tp x ,Lp y ) on S v : 

Cx = fx + Vxfy, 
Vt ~ IVVx =<Py- Vxfx- 

Solving these equations for (ip x , f y ) we find the following relations on S^: 

(1 + vl)fx = L - VtVx (1.20a) 
(1 + vl)<Pv = Vt + Vx£x ~ 7Wx- (1.20b) 

Furthermore, using (|1.20p in the boundary condition (|1.19d[) on S v we find: 



+ + + 2(ITt/I) = (L21) 

It is reassuring to note that (|1.2ip reduces to (|1.8[) when 7 = 0. 
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1.3 Multivalued Free surface 



In the above analysis we have assumed that rj(x,t) is a single valued function. 
However, in general, this assumption can be violated. In particular, it is known 
that in the case of water waves with constant vorticity, r](x, i) quickly becomes 
a multi- valued function, as observed by [T3]. In what follows, we show that it is 
conceptually straightforward to modify our analysis so that it can address the 
case of a multi- valued free surface. 

The free surface is assumed to be a one-dimensional C 2 -differentiable mani- 
fold, so features such as cusps and self-intersections are ignored. We set: 

S = {{x,y) G R 2 : x = X(X,t),y = Y(X, t)}, t > 

where A e R is the arc-length along the curve and X(X, •), Y(X, •) are assumed 
to be C 2 . The outward normal is now given by Ng = (—Y,X), where the dot 
denotes differentiation with respect to A. The kinematic boundary condition in 
I ll.l9c|) now becomes: 

(fx -IV- X t , (fly - Y t ) -N s =0, on S. (1.22) 

Again, we introduce a function £ to represent the pseudo-potential on the free 
surface, i.e. £(A, t) = cp(X(X, t),Y(X, t),t). In analogy with equations (|1.6|) . we 
note the following relations on S: 

£ = Xip x + Yip yi (1.23a) 
€t=Xt<p x +Y t <p y + <pt. (1.23b) 

Equations (|1.22p and (|1.23[) constitute a non-singular system for the functions 
{(fx, fyiPt} on the free surface S, so we can invert these relationships to express 
these functions in terms of {X, Y, r/}. After some algebra, one finds the following 
relations on S: 

^ = i*+nn ? Y + x t) -xY t) 

(X 2 + Y 2 ) V ; 

Vy = iy-x(x{ ? Y + x t) -xY t ) 

y (X 2 + Y 2 ) V ' 

t £(XX t + YY t ) + (X t Y - Y t X)(Y("fY + X t ) - XY t ) 

^ = 6 ppq^ ( L24c ) 

Using these expressions on 5, the left hand side of the Bernoulli condition in 
(|1.19d|) becomes: 

t, v lyJ iy2 , l Y ( 2XY Y t - 2X 2 X t - 2£X + 1 X 2 Y) (£ - XX t - YY t ) 2 



2(X 2 +Y 2 ) 2(X 2 +Y 2 ) 
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and the term due to surface tension is: 



XY-YX 
° {X 2 + Y 2 f/ 2 ' 

Here the coefficient of a is the intrinsic curvature of the surface S. Combining 
these expressions gives the Bernoulli condition on S in terms of the surface 
parameters (X,Y) and the potential on the free surface, £(A,i): 



Zt+gY -\X 2 t -\Y 2 



jY{2XYY t - 2X 2 X t - 2£X + ^X 2 Y) 
2{X 2 + Y 2 ) 
- XX t - YY t ) 2 XY - YX 



-a—. , -=0. (1.25) 

2(X 2 + Y 2 ) (X 2 + Y 2 ) 3 / 2 V ' 

We note that in the case X = A, Y = Y(X, t) this equation reduces to (|1.2ip . 

The above analysis shows that although it is straightforward to incorperate 
the case of a multivalued free surface, the relevant formulae become more com- 
plicated, so for convenience we will present most of our results assuming that 
T](x, t) is single valued. 



2 The Non-Local Formulation 

A novel non-local formulation in two and three dimensions for irrotational water 
waves was presented by [1] . Here we extend this formulation to rotational water 
waves with constant vorticity in two dimensions and for irrotational water waves 
in an arbitrary number of dimensions. 

The formulation by 1 is based on the existence of the so-called global relation 
(see 0). The global relation is a consequence of the following fact: suppose that 
the functions u and v are harmonic in fi C R" _1 x R. Then, 

Div {(d y u)d x v + {dyv)d x u} + d y {{d v u){d y v) - d x u ■ d x v] = (2.1) 

for each (x,y) G Q. This can be verified by expanding out the left hand side of 
(|2~Tj) to find: 

(dyv)Au + (d y u)Av, 
which vanishes in Q since both u and v are harmonic in f2. 



2.1 Irrotational n-dimensional Case 

We assume <f>{x,y,t) and r](x,t) have sufficient decay as |x| — > oo for each 
(y,t) in order for the integrals that follow to exist. This is indeed the case, 
if for example V</> € L 2 (Jl) and r\,r\t € L 2 (R™ _1 ); the justifications of these 
assumptions is beyond the scope of this present work. 
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Suppose k e R™" 1 and define v £ C°°(r>) by: 

v(x, y) = exp(ifc • x + ny), k — ±^Jk\ + ■ ■ ■ k^^. 

Then v is harmonic in Q and the following holds for (x, y) £ f2: 

Div {v(ikd y (j) + nd x (j))} + d y {v(nd y 4> — ik ■ d x <j))} = 

Integrating this expression over f2 and applying the divergence theorem gives: 

e ik-x+K V ^g^ + _ ik . g^ . ^ ^ 

+ / e ik - x+KV (ikd y (j) + Kd x <f>, Kd y (f> - ik ■ d x <f>) ■N B dx = $ (2.2) 



We have discarded the contribtions from using our assumption about the 

decay of the fields. The contribution from S v in (|2.2|) . is given by the following 
integral: 

eik . X + K y ^ ^ _ g^ g^ _ . fc _ ^ + (Q^Q^ ] fa 

From (|1.3cp the first term in the above integrand becomes %, and by (|1.6p . the 
second term becomes d x q. This gives the expression 

jk-X+Kr, ( Kr]t _ if; . g^ fa ^3) 

Now we consider the contribution from Bh in (|2.2p . We introduce Q(x,t) — 
4>(x,—h^ — h, t), which represents the potential on the bottom Bh.. In analogy 
with (|1.6p . an application of the chain rule yields: 



d x Q = d x <t>-(d y (l>)d x h. (2.4) 
Recalling that Ng = (~d x h, — 1), the contribution from Bh in (|2.2p gives 

e ik- X+K y ^ (g^ + g x(f) . g xh) _ ^ . ^ _ ( dy(j) ) dxh ) ] fa ( 2 . 5 ) 



Equation (|1.3bp implies that V(f> ■ = on Bh, so the first term in the above 
integrand vanishes. Using the result in (12. 4p we find that the expression in (|2.5p 
becomes: 

r e ik - x -^ ho+h \ik-d x Q)dx. (2.6) 



Now combining the results in ()2.3|) . (|2.6p and (|2.2p we find that equation (|2.2p 
reduces to the following equation: 

ik - x+KV {KT) t - ik ■ d x q) Ax + ( e *-x-K(ho+h) ( ik . q x q^ dx = Q 



By evaluating (|2.7p at ±|k| and adding/subtracting the resulting equations, we 
arrive at the following result: 
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Proposition 1 (Irrotational Waves in n-Dimensions) . The boundary value prob- 
lem given in (|1.3a[l - (|1.3c[l is equivalent to the pair of integro- differential equa- 
tions: 

J e lk " x [nr/t sinh(K77) — ifc ■ d x qcosh(m]) + ifc ■ d x Q cosh(K(/i + ho))] dx = (2.8) 
and 

J e lk " x [urit cosh(ft7/) — ifc • d x q w\h.{Krj) — ifc • d x Q sinh(«;(/i + ho))] dx = (2.9) 

valid for each k G R" -1 . These equations, together with the Bernoulli condition 
(|1.8[) constitute three equations for the three unknowns (t), q, Q). 



2.2 The Rotational two-dimensional Case 

Let n = 2 and also confine attention to the case in which Bh = Bq is constant. 
Using a similar approach to the rt-dimensional, irrotational case, we find the 
non-local integro-differential equation for the pseudo-potential ip(x,t) and the 
wave height n(x,t): 

e ikx±ky + ± _ ^ dx+ e ikx±ky ^ _p fa _ 

S v J Bo 

which is valid for fc € R. Invoking the boundary conditions (| 1 . 19b|) , (|1.19cj) and 
(|1.20p this expression becomes 

J e lkx±kr > [i& ± (77777, - Vt )]dx + J e ikx ^ ha iip x (x, -h ,t)dx = 0. 

Subtracting the above two expressions eliminates the Bq integral completely 
and we are left with the global relation for two dimensional water waves with 
constant vorticity 7: 

e ikx sinh (fc(jj + fc)) + i(?7t - 71%) cosh(fc(?7 + h))] dx = 0, 

where we have dropped the subscript on h. This leads us to the following 
proposition. 

Proposition 2 (Two Dimensional Water Waves with Constant Vorticity). The 
boundary value problem in (|1.19a[) - (|1.19d[) for (ip,rj) is equivalent to the follow- 
ing pair of integro-differential equations for (£, rf): 



e ikx [6s sinh (k(r) + h)) + i(r) t - 77777,) cosh(fc(?7 + h))] dx = 0, (2.10a) 

6 + & + gv + 25gggLZ%L^ - fa+^y 2 = /(„), (2.10b) 



where k G R, i) = ip(x,r](x,t),t) is the pseudo-potential evaluated on S n 
and f(rj) is defined in (|1.4|) . 
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A similar formulation can be developed if S is allowed to become multivalued. 
Indeed, if we assume S is a differentiable manifold embedded in R 2 via x = 
X(X,t) and y = Y(X,t), with \Y\ — > at oo, then a calculation similar to the 
previous result gives: 

Proposition 3 (Two Dimensional Water Waves with Constant Vorticity and 
Multivalued Free Surface). Let the free surface S be a 1- dimensional C 2 -manifold, 
and let its embedding in R 2 be parameterised by X(X,t) and Y(X,t). Then the 
boundary value problem in (|1.19a[l - (|1.19d[l for (<p,S) is equivalent to the follow- 
ing pair of integro- differential equations for (£, X, Y): 

J e [kx [i sinh (k{Y + h))+ i(XY t - YX t - jYY) cosh(fc(F + h))] dX = 0, 

(2.11a) 

jY(2XYY t - 2X 2 X t - 2£X + ~fX 2 Y) 



tt+gY- \X 2 - W 2 -r : : 

2 * 2 t 2(X 2 + Y 2 ) 

(1 - XX t - YY t ) 2 XY - YX 

ks — tj_ _ a _ ; = Q 

2{X 2 + Y 2 ) {X 2 +Y 2 f/ 2 

(2.11b) 

where k £ R, £(A,i) = ip(X(X,t),Y(X,t),t) is the pseudo-potential evaluated on 
S. 



3 A Rigorous Derivation of the Linear Limit of 
n-Dimensional Irrotational Water Waves 

Here we concentrate on the nonlinear boundary value problem described by 
(|1.3a[) - (|1.3cp and (|1.8[) in the case Bh = Bq is constant. All derivatives are 
understood in the weak sense and we work with 5'(R n_1 ) D L 2 (R™" 1 ), the 
space of tempered distributions on R" _1 . When we assert that u S iJ s (R™~ 1 ) 
it is to be understood in the sense that: 

(1 + K 2 )k(fc,-) G L 2 (R n_1 ), 

where u £ 6>'(R" _1 ) is the Fourier transform of u. The Sobolev spaces are in 
terms of the x-coordinates, and u € H B means that < oo for each 

fixed t. 

Evaluating the integro-differential equations in Proposition [T] for h = 0, i.e 
Bh constant, we find: 

/ K sinhM - i* • 8 xQ cosh( K?? ) + ifc • d x Q cosh^o)] dx = (3.1) 
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and 



J e lk ' x [nrjt cosh(nri) — ik ■ d x q sinh(K7/) — ik ■ d x Q sinh(Kh )] da; = 0. (3.2) 

Now multiplying (|3.1[) by sinh(/t/i.o) and (|3.2[) by cos!i(k/io) and adding, we find 
the following equations: 



3 [k ' x [nrj t cos!i(k(?7 + h )) - ik ■ d x q sinh(K(n + ho))] dx = 0. (3.3) 



This integro-differential equation is valid for k G R n_1 and constitutes the 
global relation for the problem in (|1.3a[) - (|1.3c|) in the particular case where Bh 
is constant. Our aim is to make suitable approximations in (|3.3p and to bound 
the relevant errors. 

The linear limit is found by assuming (77, q) and certain derivatives thereof 
are small, in an appropriate sense, and discarding terms that are smaller. It is 
convenient to work on iJ^R"" 1 ) n L^R™ -1 ) C ^'(R"" 1 ) with the following 
assumptions: 

for small e. This assumption corresponds to solutions with small energy. 

We now concentrate solely on the first term in the first integral in (|3.3p . 
since the results for the other two terms can be derived analogously. 

Lemma 1. Let fj denote the Fourier transform of n £ ^'(R™^ 1 ) n L°°(R ra_1 ) 
and suppose max{\\i]t\\ l 1 , WvWl^} < £- Then the following estimate holds: 

^. x mco S H K (v + ho)) dx= 2 

cosh(K(e + h jj 
valid for k < 0(l/e). 

Proof. We breifly outline how to prove the following basic estimate: 
ik-xVt cosh(K(?7 + h )) 



. ■ dx — rj t 
cosh [k{6 + ho)) 

First we note that the LHS can be written as: 
1 



< 2(1 



e- 2Ke ) e. 



Ak-x 



' x [rjt cosh(Ac(?7 + ho)) — rjt cosh(^(e + ^o))] dx 



(3.4) 



cosh («;(e + ho)) 

Given that ||r/||^°o < e, the following identity holds almost everywhere in R" : 



cosh (k(t]+ h ))— cosh (k(£+/i )) = e 



K(h +e) \ l e )_ 



+e" 



-n(h +e) _ 
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Using this in (I3.4[) wc find: 



' c ik-x Vt cosh(K(?7 + h )) ^ 
cosh (n(e + ho)) 



Vi 



o n(h +e) 



< 



cosh (n(h + e)) 



kl(l 



- K l"- £ l) dec 



<2(l-e-^)|MU 1 . 
Using the a priori bound H^Hl 1 < e, the result follows. ■ 

Remark 1. We can make the bound in lemma [T] sharper by assuming that 
rjt G i? s (R n_1 ) (s > 1) with || ryt 1 1 /f» < e, and then use integration by parts s 
times. This would improve our bound by some algebraic order in k. However, 
the bound is not uniform in n. The asymptotic estimate in lemma [T] is valid 
for k < 0(1 /e), which means that the estimate is valid for sufficiently large 
wavelengths. 

Remark 2. If we strengthen the condition on r\ so that n G iiP(R" -1 ) with 
s > |(n — 1), then the standard Sobolev embedding result (see for example 
ffB): 



C = (7(a). 



gives an alternative choice for the a priori bound used in the previous lemma. 
In this case it would be enough to assume ||??||h"s < e for the estimates to hold, 
but we choose to assume the bound on ||7?||l<=° since this appears more relevant 
from physical considerations. 

The second term in (13.31). i.e the term 



e k ■ d x q sinh(K(ry + ho)) dec, 



can be estimated in an entirely analagous fashion with the previous result. 
Indeed, using a similar argument and the a priori bound HfelU 1 < e wc find: 



ik ■ d x qsmh(K,(ri + ho)) 



— ik ■ d x q tanh[K/io] I dec 



< 2 1 



e- 2eK |e 



cosh(«;(e + ho)) 

Then using d x i— > ik in the Fourier integral, we arrive at the following result 



)' 

(3.5) 



Lemma 2. Let q denote the Fourier transform of q G 5'(R' i_1 ) such that 
q x G i 1 (R" _1 ). Suppose max{||77||ioo , ||<7 x ||£i} < e. Then the following estimate 
holds: 



\k-x 



ik ■ d x qsmh(K(ri + ho)) 



cosh(K(/io + e)) 



dec 



-k tanh[Kho]q + 0(e 2 ), 



valid for n < 0(1/ e). 
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The results in lemmas [T] and [2] provide us with a rigorous linear reduction of 
3.3|) valid for sufficiently long wave lengths, as summarised in the following. 



Proposition 4. Let (q,rj) satisfy the boundary value problem in (|1.3a[) - (|1.3cp 

for the case of the fiat bottom Bh = Bq. Assume that: 

max-dl^Hi,!, IMIioc, ||«x||z,i} < e- 
Then the following estimate is valid: 

fjt — Ktanh[K/io]<7 = 0(e 2 ). (3-6) 

for k < 0(1/ e). 

Proof. As a consequence of the global relation (|3.3[) , we observe the identity: 



dx), (3.7) 



ft - «t«nh[«fco]i = ( ft - / ^' ^^ d, 
1 ■ cosh(K(ft + e)) 



Ktanh[K/io]g 



ifc • d x qsmh(K(r) + ho)) 



cosh(n(h + e)) 



which holds for all k £ R™ 1 . The result follows from an application of lemmas 
[Hand EH ■ 

In what follows we will use the same contraints on (77, q) and linearise the 
Bernoulli equation in the Fourier space, so that we find an additional equation 
that couples rj and q, or rather i) and q. 

Recall that the Bernoulli condition reads: 

, 1 1» 12 , (fit + d x r) ■ d x q) 2 a I d x r) \ 
qt + ^\d x q\ +gr) = -d x ■ . (3.8) 

2(1 + \dxVn P \y/l+\d x r)\ 2 J 
It is convenient to rewrite (|3.8p as follows 

q t +gr)--Ar) + N(q,ri)=0, (3.9) 
P 

where A is the standard Laplacian on R™ -1 and N(q, rj) is defined as: 

By applying the Fourier transform to (|3.9[) it can be shown that under the 
assumption that ||f/||_ff 2 and ||<7 x ||l2 are sufficiently small, N(q,rf) is negligable 
so that the linear terms in (|3.9p are a good approximation for the dynamics. 
Indeed, using similar estimates to those in lemmas [1] and [2] (see appendix) we 
arrive at the following result. 
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Proposition 5. Let (q,r)) satisfy the Bernoulli condition Q3.8[) . Assume that: 



max {||mllL 2 , WvWh 2 , \Wx\\l 2 } < e< 
Then the following estimate is valid: 



g t +.g ? )+-K 2 ry-0(e 2 ) (3.11) 
P 



uniformly in K. 



Remark 3. The results from Propositions |4] and [5] yield the classical dispersion 
relation for linearised water waves. Indeed, differentiating (|3.6|) with respect to 
t and using (|3.11j) we find the following equation for t): 



fjtt + Kg tanh[/t/io] ( 1 H k 2 I f) = 0, 



9P 

where we have discarded the 0(e 2 ) terms. 



4 Formal Asymptotic Results for the two Di- 
mensional Rotational Case 

In this section we non-dimensionalise the equations in (|2.10[) and approach the 
problem perturbatively. Throughout this section we make the assumption that 
each of {rj, r] t ,r] x , £, £ x } are bounded and have sufficient decay at oo so that the 
results that follow remain valid. The rigorous justification of these results should 
be achieved using similar arguments to those in §3, but this is not persued here. 
We suppose £ is a typical length scale for the wavelengths and a is a typical 
amplitude of oscillation. Then we make the following substitutions: 



x^>lx, k^-, t^-j=t, £t-^-^=£, n^a-n, 7 1-* — —7. 



We introduce the dimensionless parameters (e, S) defined by: 

a h 

which are assumed to be small. In this case (|2.10a[) becomes: 

ikx {£ x sinh[£fc(e?7 + 1)] + iS[r) t - eSjr]r) x ] cosh[Sk(er] + l)]}dx = 0. (4.1) 



e 



It is straightforward to individually dominate the terms appearing in the inte- 
grand, assuming appropriate bounds on ||^x||l 2 5 Ik/ll-ff 1 an< i ll^tlU 2 - An appli- 
cation of the dominated convergence theorem allows us to expand the relevant 
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expressions as a power series in (e, S), so (|4. 1[) yields the following expression: 
f2 t n S m J e ikx A nm (k, V , T H ,£)dx = 0. 

n,m— 

Using the correspondence between k i— ► id in the Fourier integral for a finite 
number of terms (so the relevant expression is well-defined in a classical sense) 
the same equation yields the following: 

finite oo „ 

J2 € «5 m / e ik *A nm (id, r), Vt7 0dx + J2 £ n t m / e ikx A nm (k, rj, Vu da; = 0. 

n,m— n,m 

It is straightforward to bound the terms in the latter integral so the sum is 
0(e M 8 N ) for some specified M, N. Using the completeness of the Fourier trans- 
form we deduce: 

finite 

e n S m A nm (r,,r, t ,^~0, 

n,m— 

where for convenience of notation we have droppped the id dependence. The 
first few A nm can be easily computed: 

xx U 2 r Jt xx gt,xxxx 

A to>nt,Z)= o o -WntU-WZxUx ••• ■ (4 - 2) 



We note that each of the coefficients is real, and hence (|4.ip is one equation. This 
is expected, otherwise equations (|2.10l) would constitute an over-determined 
system of equations for the unknowns (t},£). 

Now we look at the equation (|2.10b[) for (rj, £). Using the non-dimensional 
parameters this equation becomes: 

6 + \<l + v 

e5"fr](2eS 2 r]r] x — 2£ x + Srj) e8 2 (i] t + er) x £ 



2(1 + ^2) 2(1+6^2) \^ 1 + £ 2 S 2 V 2 



(4.3) 



where a = a/gh 2 . Again we expand the expression in terms the dimensionless 
parameters (e, 8) to find a series of the form: 



finite 



6+ ]T e n 6 m B nm ( V ,r, t ,0~0- 



71,771=0 
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Again, the computation of the coefficients B nm straightforward 

77 -crt\xx 

_L 

2 





fit 







-ThVxZx 



(4.4) 



The above analysis yields the following systems of equations: 
Order (e°5°) : To lowest order the evolution equations are: 

£t = -v- 



(4.5a) 
(4.5b) 



Equations (|4. 5[) are a Hamiltonian system with the Hamiltonian density: 

n 00 (v,0 = h[( r i 2 + £)ax (4.6) 

and the standard symplectic structure. Indeed, (|4.5|) can be written in the form: 
d t (^=J5H 00 , J=f° 1 J). (4.7) 
Order (e 1 ^ ): To next order the equations are: 

m = -tax - e iv£,x) x , 

tt = -ri-Hx- 

This system is again Hamiltonian with the following Hamiltonian density 



(4.8a) 
(4.8b) 



ft 10 07,£) = ftoofof) + e| y ^dx. 



Indeed (14.811 can be written in the form: 



Order (e 1 ^ 1 ): To this order the evolution equations are: 

Vt = -Cxx - e(ritx)x + tSlVVx, 
St = -V- htl + e^Sx- 

The associated Hamiltonian density is: 

ftn(?7,S) =ftio(??,S) +e<5 / 77?77xSda;. 



(4.9) 



(4.10) 



(4.11a) 
(4.11b) 



(4.12) 
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Indeed, (|4.1ip can be written in the form: 



dtip =JSH n . 



(4.13) 



Order (e°<5 2 ) : At this order these exists a slight complication, because the RHS 
of the evolution equations involves r\ t : 



d t ( i) = JSHu + S 2 C 2lltxx + ^ xxxa 



(4.14) 



However, by using the expression for r\ t recursively, we can express the RHS in 
terms of (77, £) and x-derivatives thereof. One must keep track of the order of 
the relevant terms in the recursive routine. Implementing this approach gives 
to 0{8 2 ) the following equations: 



<h 1^1 =./<m,, + s- 



or\ xx 



(4.15) 



This system is again Hamiltonian with respect to the Hamiltonian density given 

by 

- Hu + \5 2 J (*vl - UL) dx. (4.16) 

Order (e 1 ^ 2 ): Again computing the relevant 0(e5 2 ) terms through the recur- 
sion process, we find the following Hamiltonian system: 



dt[')= JSHm + e&* x 



"(»7Cxx) 



The associated Hamiltonian density is: 

H21 = H20 + \e5 2 J (± 7 V - e xx rj) dx 



(4.17) 



(4.18) 



The above analysis suggests that the system is Hamitlonian at all orders, 
with the infinite chain of Hamiltonians: 



Hoo Hoi 



■ Ti-02 



H 



0:i 



H 



10 



U 



12 



H 



211 



H21 



n 



:',() 
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Conjecture. The full system is Hamiltonian with respect to the following Hamil- 
tonian: 

f ri. = CO CO 

n— m— 

where the grading is with respect to the perturbation parameters (e, 8). 

Eliminating r/ from equations (|4.15p and discarding terms of order 0(e N S M ) 
with TV + M > 2 we find the equation: 

+ eS-fWtU + Utt) + \e 2 ^ xx + <5 2 (o- - §) = 0. (4.19) 

We introduce the long length scale X and the fast time scale T by the equations 

eX = x, eT = t. 
These equations imply d x = e~ 1 dx and &t = e~ 1 dr and (|4. 19|) becomes: 

£tt — £,xx + (2£x£xt + £t£xx) 

+ 5 7 (2fr6fT + frfrr) + \&ixx + e- 2 5 2 (a - ±) Uxxx = 0. (4.20) 

Now suppose the vorticity is large, namely 7 = 0{5~ l ) and set k — jd which is 
then O(l). In addition, set 5 = e. Then (|4.20|) becomes: 

£tt — ixx + ^ix^xT + £,t£,xx 

+ k(2£t£xt + 6^tt) + f&frjr + (*-§) frxxjc = 0. (4.21) 

This equation is different to the analagous equation for irrotational waves. In- 
deed, performing calculations similar to those by pQ, one finds 

q TT - q X x + 2q x qxT + qrQxx + ^QxQxx + (o - - f ) 9xim = 0, 

where q is the velocity potential on the free surface of the irrotational fluid. 
We note that if we use in (|4.19| the change of variables 

X = x-t, r = et, 

and also set e = S 2 , then to leading order equation (|4.19| gives 

(2£ T +g + (*-i)&xx) x =0, (4.22) 

which is in agreement with a similar result derived by pQ. 

Proposition 6. Equation (|4.21[) admits two 1-parameter families of solitons 
under the following constraints: 
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1. <7 > c < 1 and K — O(l) arbitrary. 

2. a < i, c > 1 and either: 



i V2 r r .... i v/2 r r 

{%) k< \ l - or hi) k>--\ \1 -. 

c c V C c c V cr 

The cases (i) and (ii) correspond to elevated and depression solitons respectively. 

Proof. We consider travelling wave solutions in the form £(X, T) = f(z), where 
z = X — cT . Using this ansatz in (14.21|) and applying the chain rule gives the 
ordinary differential equation: 

(c 2 - 1)/" + 3c( K c 1)/'/" + §/" (f'f + / (4) = 0, (4.23) 

where the prime denotes differentiation with respect to z. We can integrate 
(14.23|) with respect to z to give: 

(c 2 - 1)/' + f («c - 1) (/') 2 + Cf T + / (3) = 0, (4.24) 

where we have assumed the derivatives of / vanish as \z\ — » oo. Now multiplying 
(|4.24p by /" and integrating once more gives 

a(f') 2 + /3(f) 3 ~ M/') 4 - (f'f = 0, (4.25) 
where we have discarded constants of integration and introduced the constants: 
, 1 — c 2 c(l — kc) 1 

a= ¥^lr p= ~^w 

For real solutions we require a 1 > 0, so we introduce set a' — a 2 . There are two 
seperate cases in which this condition holds: (i) |c| < 1 and a > | or (ii) \c\ > 1 
and ct < i. Now we introduce w(z) — f'(z), the physical variable corresponding 
to the speed of the wave propagating along S v . Rewriting (|4.25p as an ODE in 
w: 



dz\ 2 1 



Aw J w 2 (a 2 + (3w — ijlw 2 ) 
At this point it is useful to observe some obvious symmetries in (I4.26[) : 

Gi : (w,z\a,(3,n) i-> (-w,z;a,-f3,fi), 
G 2 : (w,z;a,/3,fj,) h-» (w,z;-a,f3,fj), 
G 3 : (w, z; a, /3, fi) h-» (w,z+ z ; a, j3, fi), 
G 4 : (w,z;a,j3,fi) h-> (w, -z; a, /3, /z) 



(4.26) 
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In particular, given a solution w = f(z; a, f3, fi), these symmetries generate the 
new solutions: 



(d o f)(z; a, (3, fi) = -f(z, a, -/3, fi), 
(G 2 o f)(z; a, (3, fi) = f(z, -a, (3, fi), 
(G 3 o f)(z; a, (3, fi) = f(z - z , a, f3, fi), 
(G 4 o f)(z; a, j3, fi) = f(-z, a, (3, fi). 

This will prove useful in what follows. Integrating equation (|4.26|) and invoking 
each of the above symmetries gives the solutions: 

. , , ±4a 3 exp(az) 

w(±z + z ) - 



Aa 2 fi + (3 2 =F 2a(3 exp(az) + a 2 exp(2az) 

where zo is constant. It is clear that the actions of G\ and G2 ° G4 are equiv- 
alent, so we are left with four different 1-parameter families of solutions. How- 
ever, some algebra shows that the action of G4 is equivalent to the horizontal 
translation: 



- log (4^ 



so the symmetries {Gi} only generate two independent, 1-parameter families of 
solutions: 

/ , s _ 4a 3 exp(az) 

w^z + zo) 4 a 2 /i + / j2_2 a/ 3exp(az)+a 2 exp(2az)' 1 ' 

, s -4a 3 exp(az) 

wi(z + z a) Aa 2 lM + i 3 2 + 2a peMaz)+a 2 ^ I ,{2azy { ' 

where the arrows indicate that the wave is elevated (classical soliton), or a 
depression wave. Observe that in each case the denominator is of the form: 

[aexp(±az) =p f3} 2 + 4a 2 fi. 

It is clear that in case (i) |c| < 1 and a > |, this expression is positive for all z, 
since fi > 0. So for \c\ < 1 and a > |, equation (|4.21| admits 6oi/i depression 
solitons and elevated solitons. 

In case (ii) we see that fi < 0, so equations (|4.27[) do not give soliton solutions 
when the denominator is allowed to vanish. For w-\ to give a soliton solution, 
we require (3 < and Aa 2 fi + (3 2 > 0. Using the definitions of a, /?, fi, we see 
this is equivalent to: 

c(l - kc) > 0, 
c 2 (1-kc) 2 + 2(1-c 2 ) >0. 
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Figure 1: Depression and elevated solitons corresponding to (c, re, a) = 
(0.95,0.5,0.4). 

We take c > 1 and so the first inequality gives rec < 1, and the second gives: 

«<--^ (4.28) 
c c V c 

Similarly, for uij to give a soliton solution we require (3 > and 4a 2 /x + 1 > 0. 
In analogy with (|4.28|) we find the condtion: 

+ — \A-4 (4-29) 
c c V cr 

This completes the proof of the result in Proposition [6] ■ 

We note that in the absense of vorticity (i.e re = 0), neither of the conditions 
(|4.28[) and (|4.29p can be satisfied, so the equation (|4.2ip admits no classical 
soliton solutions when a < |. This is in agreement with the rigorous non- 
existence results of [T5] . 

It should be noted that the PDE in (I4.21|) is not well-posed in general when 
a < ^ . Considering plane wave perturbations with small amplitude leads to the 
dispersion relation: 

u 2 = k 2 + (a - |) k\ 

In the case a < ^ we see that the amplitude of the perturbations will grow expo- 
nentially in time. This indicates that solitons derived in case (6) of Proposition 
[6]will be highly unstable. Nevertheless, it is interesting to look at the behaviour 
of these solitons which are completely ruled out in the absense of vorticity. It 
is convenient to perform a translation in z so that the solutions are symmetric 
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Figure 2: Depression and elevated solitons in the case a < | for K = § and 
c> 1. 



about z = 0. Some algebra gives: 

w T (z) = 



1 + T2 cosh(az) 



1 — cosh(az) 
where we have introduced the two constants: 

2a 2 „ 1 



Ti^- — , T 2 = -y/4»a?+lP. 

In fact, if we don't specify which side of the branch we are on in the (-) 1//2 term 
in T2, these two solutions can be encorperated into one single expression: 

W(z) = Fl - . , (4.30) 

w 1 + T 2 cosh(az) v ; 

which then describes elevated solitons when wc take the positive square root and 
describes depression solitons when we take the negative square root. In this case 
the constants \Ti, T2} are determined by the points on a 1-dimcnsional manifold 

r : 

r = {(r 1; r 2 ) e R 2 : r 2 - (4) r 2 = 1} (4.3i) 



In summary, the soliton solutions of equation (|4.2ip are given by £(X, T) = 
W(X - cT) where: 

w ( g ) = TTr^-Ti V ( r i- r 2)er. (4.32) 
1 + I 2 cosh(az) 

There is now a clear geometrical interpretation of the different types of soliton 
solution to (|4.21|) in terms of the curve T. The topology of T is completely 
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determined by the sign of /i/a 2 . Using the definitions of a and fi, one sees this 
is equivalent to the sign of a — | , or the sign of c 2 — 1 . 

1. (a > |). In this case T is hyperbola, i.e. two disconnected copies of 
R. There is a bifurcation on each of these disconnected componets at 
r 2 = 0, which corresponds to constant solution W = T%. Removing 
these two points we are left with four disconnected components. The two 
components with > correspond to elevated solitons, whereas the 
two components with Til^ < correspond to depression solitons. 

2. (cj < i) . In this case T ~ S 1 . In this case there are four bifurcations. As in 
the previous case, two of these bifurcations are at T 2 = and correspond to 
the constant solution W = Y\. However, at the bifurcations points I\ = 0, 
the solution becomes unphysical and has finite time blow-up. Removing 
the points of bifurcation, we are left with four disconnected components. 
The two components in r 2 > correspond to elevated solitons (ri > 
0) and depression solitons (Ti < 0). The remaining two disconnected 
components correspond to unphysical solutions which blow up at az = 
arccosh(|r 2 | _1 ). 

5 Conclusions 

We have presented the following generalisations of the results of pQ: (a) We 
have extended the formalism to the case of n-dimensions, n > 1; (b) we have 
considered the case of constant vorticity; and (c) we have incorperated the effect 
of a multi- valued free surface. 

We recall that [7] introduced an elegant Dirichlet to Neumann operator G(rj) 
associated with the velocity potential and also obtained a series for the operator 
G(r)), valid for small 77. The integral equations present in [T] and here, can be 
considered as the summation of the above series, i.e. the series of [7] is the 
Neumann series of the integral equations derived by [T] and here. 

It appears that the new formulation provides an alternative, perhaps simpler 
approach, for (a) the numerical investigation of water waves; (b) the derivation 
of various asymptotic limits; (c) the rigorous analysis of water waves. 

Regarding (a), we recall that two-dimensional lumps were computed by [1] 
in the case of sufficient surface tension. 

Regarding (b) we recall that various asymptotic equations, including the 
Boussinesq, Benney-Luke and the nonlinear Schrodinger equations, were de- 
rived by pp. Similarly, an appropriate Boussinesq type equation in the case of 
constant vorticity is found in §4. We observe that these equations have been 
derived by several other authors, however, it appears that the new formulation 
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provides a straightforward way of deriving these equations. For example, re- 
garding our results in §4, we note that [5] uses a perturbative approach which 
involves solving many PDEs arising from the consistency of the perturbation 
expansion with the boundary conditions. To solve these PDEs, the author con- 
fines attention to seperable solutions. The entire paper is devoted to computing 
the first few coefficients of the perturbation expansion. The author comments 
that this method is unweildy at anything beyond second order. Our paper gives 
a much more direct approach, and can easily be extended to higher orders. It 
is not immediately obvious as to whether his equations are exactly the same as 
ours, because our soliton equation is for the wave speed, whilst his are for the 
wave height. 

Regarding (c) we have shown in §3 that standard PDE techniques can be 
used for the rigorous analysis of water waves, at least in the linear limit. The 
new formulation suggests a rigorous methodology which differs drastically from 
that employed in the important works of [H O [TB] ; the extension of the results 
of §3 to the nonlinear problem is a work in progress. 

Finally, we note that in the classical works of [I] , the author focuses entirely 
on the amplitude rj{x, t) and does not include the effect of surface tension. In 
contrast, the solitons discussed in this work concern the wave speed, £(x,i), and 
our formulation also includes the effect of surface tension. Our analysis, partic- 
ularly the results in Proposition [6l outline the importance of the relationship 
between the surface tension and the vorticity for the existence of solitary waves. 
Furthermore, our formalism can easily be extended to the three dimensional 
case. The main advantages of the new approach presented by [1] and in this 
paper, are a consequence of the explicit nature of the equations in Propositions 
III land [3] 
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6 Appendix 



Here we prove the result in Proposition [5j It suffices to prove the estimate: 

\\N(ri,q)\\ Ll = 0(e 2 ) 
given that max{||ry|| ff 2, < e. Recall that: 



N(q,r ] ) d ?t\d xq \ 2 + -d x 



P \^l + \d x r,\* 
The following estimate is clear from the definition 



dxV 



(r]t + d x i] ■ d x q) 2 
2(l + \d xV \ 2 ) 



\\e ik - x N( V ,q)\\ L1 <l\\d x q\\l 2 +^ 



\/l+\d*v\'- 



- d x r] 



) 




nt+dxTl-dxq 




L 1 


x/i+is.')! 2 



L 2 



The first term is clearly 0(e 2 ) since < e. For the third term, we use the 

estimate: 



< IMU' + 



\d x q\ 



\d*v\ 



2 



L- 



<(IMU» + RgM (6-1) 
which gives the required 0(e 2 ) bound. For the second term, note that: 



d x 



d x r) 



d x r] 



< 



L 1 



1 A?7 



+E 

L 1 i,3 



(i + |d^| 2 ) 3 / 2 



L 1 



where r\i = drj/dxi. Rewriting the second term, we observe the following esti- 
mate: 



E 

hi 



ViVjVij 



(1 + \d xV \ 2 ) 3 ^ 



L 1 



= E 
^E 

i,3 



(i + E fe % 2 ) 3/2 



VijVj 



L 1 



(1 + ^)3/2 

^EiTS^II^ 

^EiTf II%'IU 2 IMIl=> 



L 1 



(6.2) 



where we applied the Cauchy-Schwarz inequality. Similarly, we find: 
1 



-1 At) 



< i||fl^|| ia ||A»j|| La . 



(6.3) 



L 1 
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Combining (|6.2p and (|6.3p we find: 

< Nil- (6-4) 



L 1 

From the estimates in (|6.1[) and (|6.4[) it is follows that \\N(rj,q)\\ L i — 0(e 2 ). 
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